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Abstract 

In this paper, we discuss a possibility of studying properties of dark energy in long baseline neutrino 
oscillation experiments. We consider two types of models of neutrino dark energy. For one type 
of models the scalar field is taken to be quintessence-like and for the other phantom-like. In these 
models the scalar fields couple to the neutrinos to give rise to a spatially varying neutrino masses. 
We will show that the two types of models predict different behaviors of the spatial variation of the 
neutrino masses inside Earth and consequently result in different signals in long baseline neutrino 
oscillation experiments. 


I. INTRODUCTION 

There are growing evidences from various cosmic observations, including type la supernova (SNIa) |lj, cosmic 
microwave background (CMB) Q, large scale structures (LSS) (3,0, and so on, that support for a spatially fiat and 
accelerating universe at the present epoch. In the context of Friedmann-Robertson-Walker cosmology, this acceleration 
is attributed to the so-called dark energy Q ■ The simplest candidate for the dark energy seems to be a remnant small 
cosmological constant. However, many physicists are attracted to the idea that the dark energy is due to a dynamical 
component in the evolution of the universe, such as the quintessence i 0 S1 US , the K-essence 0 M 0 , the 
phantom 0, or the quintom T7, 18, [20]. 

Recently there have been a lot of works 0 e m m 0 n mm m ns m 0 m m m m m on 0 , 
Ell It 3| which study possible connections between neutrinos and the dark energy, generally referred to as the neutrino 
dark energy. One of the predictions of the class of models of neutrino dark energy is that neutrino masses are not 
constant, but can vary as a function of space and time. This general prediction can be tested with Short Gamma 
Ray Burst l29(, CMB and LSS (33| . and much more interestingly and directly in neutrino oscillation experiments 
124 125L H4144 l45l l46t I 47 II . In this paper we make a concrete study of the possibility of probing the property of dark 
energy and differentiating its dynamic origin in very long baseline neutrino oscillations. 

In general for the models of neutrino dark energy, the Lagrangian is given by 

£ = + 4 + £ in t , (1) 

where C^ M is the Lagrangian of the standard model describing the physics of the left-handed neutrinos, C$ is for 
the dynamical dark energy scalar <j> such as quintessence or phantom, and £,; nt describes the sector that mediates the 
interaction between the dark energy scalar and neutrinos, and gives rise to variations of the neutrino masses. 

At energy much below the electroweak scale, the relevant Lagrangian for the neutrino dark energy can be written 
as 
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where £„ is the kinetic term of neutrinos, c is a coefficient which takes the value of 1 for a Dirac neutrino and 1/2 for 
a Majorana neutrino, rrij((f >) is the scalar field dependent mass of the j-th neutrino that characterizes the interaction 
between the neutrino and the dark energy scalar. 

The authors of have used the recently released SNIa data to constrain the coupling of the scalar (j) to neutrinos 
and the property of the dark energy scalar. They found that the model with a phantom scalar is mildly favored. 
However, the data do not rule out the possibility of the quintessence scalar coupled to neutrinos. In this paper we 
will show that these two models predict different spatial variation of neutrino masses inside Earth and consequently 
result in different signals in very long baseline neutrino oscillations. 

This paper is organized as follows: in section II we present our mechanism for the neutrino mass variation; in 
section III we study quantitatively the mass-varying effect in the long baseline neutrino oscillations; Section IV is a 
brief summary. 


II. MECHANISM FOR VARIATIONS OF NEUTRINO MASSES 

In the standard model of particle physics, a typical term for the neutrino masses can be described by a lepton 
violating dimension-5 operator 


-£v = jl L l L HH + h.c., (3) 

where / is a scale of new physics beyond the standard model which generates the B — L violations, l and H are the 
lepton and Higgs doublet, respectively. Here we neglect the lepton generation symbol. When the Higgs field gets a 
vacuum expectation value, ( H) = v ~ 174GeV, the left-handed neutrino receives a Majorana mass ~ y-. In Ref. [23j 
the authors proposed a model where the dark energy scalar <f> couples to the dim-5 operator. In this model neutrino 
masses vary along with the evolution of the universe and the neutrino mass limit imposed by baryogenesis is modified. 

The dimension-5 operator above is not renormalizable. It can be generated from physics beyond the standard model 
which involves very heavy particles interacting with standard model particles. At low energies the heavy particles can 
be integrated out and thus resulting in effective, nonrenormalizable terms. For example, in the model of the minimal 
see-saw mechanism, we have the neutrino mass term, 

~£ = Y1 h ijiLiHv R j + 2 5Z ' U o"/o"/0 + h - c -> ( 4 ) 

ij ij 

where My are the Majorana mass matrix elements of the right-handed neutrinos and /ly the Yukawa couplings. The 
Dirac masses of the neutrinos are given by muy = hjj'u. Now integrating out the heavy right-handed neutrinos VRj, 
one will generate a dim-5 operator as stated above. As pointed out in Ref. [2fj, there are various possibilities to have 
the light neutrino masses varied, such as by coupling the quintessence field to either the Dirac masses or the Majorana 
masses of the right-handed neutrinos, or to both. 

In this paper we consider the case where the variation of the neutrino masses is caused by a coupling of the dark 
energy scalar </> to right-handed neutrinos. With the Majorana masses of right-handed neutrinos varying, My becomes 
a function of the dark energy scalar field (j), My = My (</>). Furthermore, we assume a linear relationship between the 
Majorana mass and (j). Then the relevant Lagrangian can be written as 

-£ = hijl Li Hv Rj + -gij&VRiVRj + h.c. =F + V (/>), (5) 

where py are dimensionless constants and V (</>) is the potential for </>. The upper minus sign in is for the case of 
quintessence while the lower plus sign for the phantom. This convention will be used throughout this paper. Now the 
Majorana mass matrix elements of the right-handed neutrinos can be written as 

Mij = gij(f> ( 6 ) 

and consequently via the seesaw mechanism we obtain the masses of the light neutrinos v = vr + , 

1 


m u oc 


( 7 ) 
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Similar to the study on the mass varying neutrinos in |25| . here we introduce also a coupling between (p and the 
baryon matter with the effective potential for <p at low energies given by 2H 23 


V eff (cp) = V(cp) + ^ Pi e 0i ^ m ^, 


( 8 ) 


where /3i is a dimensionless constant 001, pi denotes the energy density of the i-th matter field, and m p i is the 
reduced Planck mass. The dark energy scalar cp shall change its value in space 2H 23 according to the equation of 
motion 


= ±{V,t + V 

^ rn p i 


(9) 


In the following we will calculate the evolution of the dark energy field and the corresponding variation of neutrino 
masses in the Earth for both the quintessence and the phantom cases. 

The density profile of baryon in the Earth is taken as the widely adopted PREM model 23, in which Earth is 
taken to be spherically symmetric. The atmosphere is treated as a homogenous layer of 10 km in thickness with a 
constant density p a tm — 10 -3 g/cm 3 . Defining x = r/R® and R a tm = R® + 10 km with R® = 6371 km the Earth’s 
radius, the baryon density can be expressed as 


P i(r) = at + bix + dx 2 + diX 3 for r»+i < r < n 


with 


(a u a 2 , ..., an) = ( p a tm , 1-02,2.6,2.9,2.691, 7.1089,11.2494,5.3197, 7.9565,12.5815,13.0885), 

(&n fo 2 , -, fon) = (0,0,0,0,0.6924, -3.8045, -8.0298, -1.4836, -6.4761, -1.2638,0), 
(ci,c 2 ,...,cn) = (0,0,0,0,0,0,0,0,5.5283,-3.6426,-8.8381), 

(di, d 2 ,..., dn) = (0,0,0,0,0,0,0,0,-3.0807,-5.5281,0) 

in units of g/cm 3 , and 

(n, r 2 , ..., n 2 ) = (Ratm, R®, 6368,6356,6346.6,6151,5971,5771,5701,3480,1221.5,0) 
in units of km. We also assume a homogeneous baryon background outside the atmosphere with the density 

p§ ~ 4% x p c ~ 1.4 x 10 _29 g/cm 3 , 

where p c — 4.1 x 10 -47 GeV 4 is the critical energy density of the universe at the present epoch. 

With these assumptions, Eq. © can be simplified as 

As 


d 2 cj) 2 d(f> dP 

dr r dr L dip 


m p i 


-pBe 


\B<fr/m p i 


} 


with the boundary conditions 


( 10 ) 


(ID 

( 12 ) 

(13) 

(14) 


(15) 

(16) 


= 4>u for r = r c , 


(17) 


(18) 


dd 

— = 0 for r = 0. 
dr 


(19) 


Here r c denotes the interface between the static solution and the cosmological one. For r > r c , we expect <p to become 
a constant: cj> = cpu and (pu is the value of the dark energy scalar field on cosmological scales at the present epoch. 
For simplicity we take r c ~ R a tm which is consistent with the assumption in Eq. m that the baryon background 
becomes very thin and homogeneous for r > R a tm■ It should be noted that we require A b < 10 -4 to satisfy the 
equivalence principle constraints 23 01 ■ 

In addition, we take 


T rr i\ 




/ on \ 
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as an example 15.'if and then write Eq. El as 


= ±{—— Voe-W™* 1 + —p B e 
dr z r dr m p i m p i 


For this given potential, the value of (j)u is 


m p i 2 Vo 

<t>u = —f- ln[- 


P (1 ~w)P<i> 


(21) 


( 22 ) 


with w = and + V for the quintessence case, while w = \+ v an d P<t> = — \4> 2 + V for the 

phantom. Here w ~ —1 is the equation of the state and p<f, ~ 73% x p c ~ 3.0 x 10 -47 GeV 4 is the energy density of 
the dark energy at the present time 1 . 

Assuming (f>/m p i <C 1 for 0 < r < R a tmi which we will show in the later numerical results is reasonable, we can 
simplify Eq. <ED as 


d?(j> 2 d(j) P 1 X B x 

~T + ~T = ±1-~ +- P B s 

dr r dr m p i 2 m p i 


(23) 


and obtain the solution 


<t>u 


9i-%± {-^K 1 - w )P* R lf + £nX R l(f x2 + ft* 3 + tF + M* 5 )} for r ‘+s <r<r., 


2 x 2 A b o2 (ai 2 I bi 3 i Ci „4 d* 5\ 


for r > R n 


(24) 


with 


fi—\ f. I As. c»2 / n.+i-qj 3 I fcj+1-fe.- 4 , cj+j-Ci 5 . rfj+i-dj 6 \ f nr i < 1 f) 

A+ 1 1)1 it ©l 3 ^i+l + 4 x i+l' 5 x i+l + 6 ^i+ll 101 1 — 1U 


for i = 11 


m p i 


(25) 


and 


= J <fo + Ir ± {^ 1(1 “ u FF©ir - + iH + fFi + §*1)} for * = 1 

91 =F f + for i > 2. 


(26) 


Here we have adopted the definition of xp Xi = r,/i2®. 

In the numerical calculation, we take P = 1 and <fru ~ 10 _15 m p j by choosing Vo to satisfy the cosmological 
observations, p ^ ~ 3.0 x 10“ 47 GeV 4 and w ~ —1. In Fig. |T| we plot () as a function of r with A B = 10 5 and 
10~ 6 for the quintessence and phantom cases, respectively. Our results show that in the two cases the variations of 
the dark energy field can be sizable inside the Earth. Accordingly, as shown in Fig. |2j the neutrino masses could 
have a significant variation inside the Earth. Meanwhile we notice that the consequences of the quintessence and the 
phantom are different due to the opposite behaviors of the dark energy fields. 


III. MASS-VARYING EFFECT IN LONG BASELINE NEUTRINO OSCILLATIONS 


Now we discuss the mass-varying effect induced by the evolution of the dark energy in neutrino oscillations. The 
neutrino propagation in matter is governed by the Shrodinger equation: 


i 


d 

dt 





0 0 0 
0 A rri2i 0 
0 0 AjtF 


(N e (x) 0 

U^+V2G f \ 0 0 

\ 0 0 



(27) 


1 Here we do not expect the dark energy model to preserve the tracking behavior HE1 and in this sense, the dark energy scalar field is 
essentially like the cosmological constant on the largest scales. This is also consistent with the static condition adopted in Eq. 
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FIG. 1: The evolution of the dark energy held <j> in the baryon matter background with A b = 10 -5 and 10~ 6 for the quintessence 
and the phantom cases, respectively. It is shown that the assumption of 4>/m p i 1 can be fulfilled in the whole space. P and 
Q are for the quintessence and phantom cases, respectively. 


Here U is the usual 3-flavour vacuum mixing matrix and \/2 GfN £ (x) is the MSW term (3IH]- We have also adopted 
the definition that <p s , A mfj are the values of the dark energy field and the neutrino mass squared differences on the 
Earth surface, respectively. In the following numerical estimate, we will take Am^i ~ 7.9 x 10 -5 eV 2 , tan 2 6*12 ~ 0.4 
from KamLAND [5(1, ~ 2.8 x 10 -3 eV 2 , sin 2 2023 — 1.0 from K2K j^j, sin 2 2013 <0.1 from Chooz [ 53 ] which 

we will take to be zero for simplicity, and a zero Dirac CP phase in Eq. El). 

In Fig. El by taking the longest baseline L = 2 Rq which appears in the atmospheric neutrino oscillations and 
L = 732 km which is Fermilab to Soudan or CERN to Gran Sasso (right panel), we plot the survival probabilities of 
with A b = 10 -5 , 10 -6 and 0 for the case of the quintessence and the phantom, respectively. For As = 0, the spatial 
variation of the dark energy is negligible and the neutrino oscillation is identical to the case of decoupling between the 
neutrino mass and dark energy. For A b = 10 -5 , it is clear that the survival probabilities differ significantly for the 
cases of the quintessence, phantom, and decoupling of neutrino and dark energy for L = 2 R$. But for L = 732 km 
the different cases can not be distinguished. 

We also consider the very long baseline L = 7332 km which is the distance from Fermilab to Gran Sasso underground 
laboratory in Italy j5fj and L = 9400 km which is the distance from Fermilab to Beijing. As shown in Fig. QJ the 
survival probabilities P ^ with A b = 10~ 5 is sensitive enough to distinguish the cases of the quintessence, phantom, 
and decoupling. The most sensitive measurement is perhaps at the second zero of the survival probability. For 
L = 9400 km, the decoupling case has the zero at about muon neutrino energy of 20 GeV, the quintessence above 20 
GeV, and the phatom below 20 GeV. The separation in energy is sufficiently large that it should make the distinction 
of the three case clean. 
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FIG. 2: The variation of the neutrino masses in the baryon background with As = 10 -5 and 1CP 6 for the quintessence and 
the phantom cases, respectively. Here m is the masses of the left-handed Majorana neutrinos and m a denotes its value on the 
Earth surface. P and Q are for the quintessence and phantom cases, respectively. 




E (GeV) E (GeV) 


FIG. 3: The survival probabilities P of in the long baseline L = 2I?® and L = 732 km with As = 10~ 5 , 10 -6 and 0. 
P and Q are for the quintessence and phantom cases, respectively. Note that A /3 = 0 is identical to the case of decoupling 
between the neutrino mass and dark energy. In the calculations, we have used ~ 7.9 x 10 -5 eV 2 , Am^ ~ 2.8 x 10 -3 eV 2 , 

tan 2 0i2 — 0.4, sin 2 2023 — 1-0, sin 2 2013 ~ 0, and a zero CP phase. The left panel if for L = 27?® and the right panel L = 732 km 






















7 




FIG. 4: The survival probabilities of in the very long baseline L = 7332 km and L = 9400 km with A b = 10 -5 , 1CU 6 , 
and 0. P and Q are for the quintessence and phantom cases, respectively. Note that \p = 0 is identical to the case of decoupling 
between the neutrino mass and dark energy. In the calculations, we have used AmJi — 7.9 x 10 -5 eV 2 , Am^ ~ 2.8 x 10 -3 eV 2 , 
tan 2 0i2 — 0.4, sin 2 2023 — 1.0, sin 2 2013 ~ 0, and a zero CP phase. The left panel is for L = 7332 km and the right panel for 
L = 9400 km. 


IV. SUMMARY 

In this paper, we discuss a possibility of studying the dark energy property with long baseline neutrino oscillation 
experiments. We consider two types of models of neutrino dark energy where for one model the scalar is taken to 
be quintessence-like and for another model phantom-like. These scalars couple to the neutrinos which give rise to a 
variation of the neutrino masses. We take a specific scalar dark energy potential with the Earth baryon background 
and then calculate the spatial variation of the dark energy field for the case of the quintessence and the phantom, 
respectively. We find the corresponding evolution behaviors of the neutrino masses inside the Earth could be sig¬ 
nificantly different in the two cases and hence the property of the dark energy may be probed in the long baseline 
neutrino oscillations. 
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